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SUMMARY 


A reevaluation of existing flight data obtained for a supersonic 
canard missile conf iguration of the boost-glide type has been made to 
determine derivatives not previously evaluated. These derivatives, 
together with those previously evaluated and published, have been 
utilized to determine some typical airframe frequency responses based 
on the three-degree-of -freedom longitudinal equations of motion. 

For constant flight conditions, it is shown that besides the usual 
drag-coefficient variation, the velocity derivatives also vary with 
angle of attack or Mach number. The effect of angle-of -attack variation 
on the missile frequency responses which include the velocity derivatives 
in their solution is appreciable only in the low-frequency region of 
operation (below 10 radians/sec). When the velocity derivatives are 
neglected, this low-frequency variation with angle of attack is different. 
Some of these differences are pointed out in the results to indicate the 
significance of including the velocity derivatives in the solution for 
the transfer functions. 


INTRODUCTION 


Increased interest in the phugoid motion has been expressed from 
various sources, especially in connection with automatic airspeed or 
altitude control systems. The present study of the three-degree-of- 
freedom motion actually arose out of an investigation of an automatic 
altitude control system wherein the altitude- sensing instruments operate 
at a frequency considerably below the high-frequency airframe mode. In 
this situation, the influence of the low-frequency airframe dynamics on 
the overall stability of the system may prove to be important . 


^-Supersedes recently declassified NACA Research Memorandum L54C02 
by Ernest C. Seaberg, 1954. 
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The main purpose of the present paper is to evaluate the longitudi- 
nal transfer functions of a supersonic canard missile by means of a 
reevaluation of the data obtained from previous flight tests of the mis- 
sile in order to include the variation of drag coefficient with angle of 
attack and the variation of force and moment coefficients with velocity 
in the solution of the three-degree-of -freedom equations of motion. The 
variation of these additional derivatives with Mach number is presented 
herein for several values of angle of attack. The longitudinal transfer 
functions in the form of frequency-response plots for several flight 
conditions are also given. 

In reference 1, the effect of a fourth degree of freedom (namely, 
variation of height) on the longitudinal motion of full-scale aircraft 
is discussed. Although this additional degree of freedom may influence 
the longitudinal motion of missiles to some extent, the study of this 
effect was not included in the present investigation. The present inves- 
tigation is mainly concerned with the additional effects of variation of 
force and moment coefficients with velocity, as compared with previous 
three -degree -of -freedom studies (ref. 2, for example) wherein these 
effects were not considered. 


SYMBOLS 

0 pitch -attitude angle, radians 

a angle of attack, radians unless otherwise indicated 

P angle of sideslip, radians 

7 flight -path angle, radians 

£ critical damping ratio 

6 canard elevator deflection, radians 

m mass, 5*05 slugs 

W weight, 163 pounds 

S wing area, 4.1 sq ft 

Y stability axis which passes through center of gravity and 

is perpendicular to vertical plane of symmetry 

Iy moment of inertia about Y-axis, 31-3 slug-ft^ 


NACA TN 4192 


3 


c 


X 


sm 
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AV 

V 
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t 

D 

G(D) 

|g(d)| 

db 

0 


c L 


C D 

^m 


mean aerodynamic chord, 1*776 ft 
static margin, ft 

Mach number 
velocity, ft /sec 


dynamic pressure, lb/sq ft 
time, sec 

differential operator, d/dt 

transfer function expressed as a linear function of D 
amplitude or magnitude of G(D) 
decibels, 20 logqo| G ( D )| 

phase angle of G(D), deg 
angular frequency, radians/sec 
undamped natural frequency, radians/sec 
lift coefficient , Lift /qS 
drag coefficient, Drag/qS 


pitching-moment coefficient, 


P it ching moment 


Ci* - 

S> * Xjdis 




M m 
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S 

§ 

C D<x 

= OC-^y^a 

C D U 

= dC D ^du 

°Lu 

= dCj^du 

c m u 

= ^ c m l bu 


Subscripts: 


phugoid 

equilibrium value 

Dot over a symbol denotes derivative with respect to time. 


METHODS AND APPARATUS 
Equations of Motion 


The three-degree-of -freedom longitudinal equations of motion, 
assuming small disturbances from steady symmetric horizontal flight, 

(iL cos 7o )e + (c Da - JjL cos 7o )a + + 2C D ju = 0 

(§ d - & sin 7 °) e - (§ d + ci u - & sin 7 °) a - K + 2cl ) u = c 

- f/v) e " ( Cma + i/^ D ) a " = Cm&& 


are: 
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Various derivations and forms of these equations can be found elsewhere. 
(See, for example, refs. 2 to J.) In these references various notations 
are used; in some more terms are included than in the present report, and 
in others terms which are included herein are neglected^ The equations 
presented here, however, satisfactorily define the motion of the missile 
under consideration with the assumption of small disturbances. 'The nota- 
tion used is based on the stability system of axis and parallels the more 
widely used two-degree -of -freedom equations. (See ref. 8.) 


Flight -Data Presentation 

Airframe.- The airframe on which this analysis is based is an all- 
met al~research model of the canard missile type. A photograph and a 
sketch of the actual model used for previous flight testing is shown in 
figure 1. References 9 and 10 give more complete descriptions of this 
model and present the flight data in the form of aerodynamic stability 
derivatives used in previous two-degree-of -freedom analyses. In refer- 
ence 8, values of these derivatives are tabulated for typical flight 
conditions. The additional derivatives necessary for the three-degree- 
of -freedom study presented herein were obtained from a reevaluation of 
the flight data presented in references 9 and 10 and from the results of 
other flight testing not previously published. 

Cn •- The derivative C n was evaluated from the results of flight 
•Mx v a 

tests not previously reported, wherein the model was roll stabilized and 
pulsed in pitch and yaw, and the values of Cp based on total drag along 

I 2 2 

the velocity vector were plotted against total angle of attack ^ a + 3 . 
Sufficient curves were obtained to cover a Mach number range of approxi- 
mately 1 to 1.7. The slopes of tangents to curves obtained in the fore- 
going manner were then measured to determine the Cp a data presented 

herein. Since a symmetrical cruciform missile is being dealt with, a 
value of C-n based on total drag and total angle of attack is equiva- 

U CL 

lent to one based on measurements in the pitch plane alone and is valid 
for use in the longitudinal equations of motion. 

C Du , CLu , and Cm u .- In this case the variations of the coeffi- 
cients C L , and with velocity were first obtained. Sufficient 

data were available to cover a velocity range of approximately 1,000 
to 2,000 feet per second, assuming standard sea-level conditions. The 
angle-of -attack range covered was as high as 10° or 12°. The slopes of 
tangents to these curves were then measured to obtain the nondimens ional 
derivatives such as Cp u = 8Cp/8u which are presented herein as functions 

of Mach number for constant values of angle of attack. 
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Frequency Responses 

A solution of the equations of motion to obtain the transfer func- 
tions can only be obtained by assuming an equilibrium Mach number and 
then selecting values of Cp^, Cp u , Cl^* an< ^ ^u ^ or 

of -attack value under consideration. Solutions based on the foregoing 
assumptions are presented in the results as frequency-response plots of 
the transfer functions. A particular transfer function, for instance 
0 / 5 , is obtained by expansion of its determinant to yield the function 



The frequency response of this function is then plotted through the use 
of templates as described in chapter 8 of reference 11. 


RESULTS AND DISCUSSION 
Derivatives 


Cp .- Plots of Cp against Mach number for constant values of a 

are shown in figure 2. For a constant Mach number, this derivative exhib 
its a decided increase with increasing a and for the lower values of a 
Cp a decreases to a steady value with increasing M. The decrease with 

increasing M is more apparent at the higher values of a for the Mach 
number range shown. 

CDu, Ci^, and C mu .- The velocity derivatives Cp u , Ci^, and Cm u 

are plotted against M for constant values of a in figure 3* The 
trends exhibited by Cp u and Cl^ (figs. 3(a) and (b)) are somewhat 

similar in that at subsonic Mach numbers the curves show positive values 
which decrease to zero at approximately M = 1.03 and then increase 
negatively to maximum values between M = 1.2 and M = 1.35* From this 
point, a more gradual rise back toward zero is noted for Ct than for 

Cp . 

u 


Curves of C mu against M for constant values of a are shown in 
figure 3(c). From these curves it is seen that C^ is negative at sub- 
sonic Mach numbers, and is zero for all values of a at approximately 
M = 1.03. The values of C^ then increase to maximum positive values 
at M = 1.55 approximately, after which a decrease with increasing M 
is shown. The uniform spread in the curves of this figure indicates that 
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the variation of this derivative with a for a constant Mach number is 
fairly linear. 

Cp and Cp.- The values of Cp and Cl used in computing the 

transfer functions presented herein are shown in table I. Except for 
the zero-lift drag coefficient, the Cp values have not been previously 

published. The Cl values are presented for completeness although they 
are indirectly available in published form elsewhere. 


Frequency Responses 

Typical missile frequency responses obtained from a solution of the 
three-degree-of -freedom equations of motion are shown in figures 4 and 5. 
These responses are based on standard sea-level conditions at M = 1.6 
with a missile static margin of 0.294c. The responses 9/6, a/6, 

and u/5 were obtained directly from a solution of the equations of 
motion as given herein previously, whereas the substitution y - 0 - a 
was employed to obtain the y/b responses (figs. 4(d) and 5(h)). For 
each of the variables, the amplitude and phase response obtained for 
trim angle-of -attack values of 0°, 4°, and 10° are shown. This implies 
that the derivatives Cp^ , Cp u , Cj^, and C mu and the values of Cp 

and Cl correspond to the trim angle-of-attack value under consideration 
for each solution, despite the variations with a shown in figures 2 
and 5 and table I. 

Effect of including velocity derivatives .- The frequency responses 
obtained at M = 1.6 including the effect of velocity derivatives are 
shown in figure 4. Phase shifts as the frequency approaches zero and 
changes in static sensitivity are shown for the trim angle-of-attack 
values investigated. In figures 4(a) and 4(d), for example, as the fre- 
quency decreases, the phase curves for 9/6 and y/b are shown to tend 
toward -90°, 0°, and l80° as the trim angle of attack is changed from 
0° to 4° to 10°, respectively. At a trim a of 0° the solution is simi- 
lar to that obtained for two degrees of freedom with an integration 
effectively present because of the zero lift force. At a = 4° the 
characteristic equation (denominator of transfer function) is of the 
biquadratic form 

(d 2 + 2£ p o^D + o^ 2 j(l) 2 + 2^D + a* 2 ) 

with a positive static sensitivity yielding a zero initial phase relation. 
For a = 10°, however, the numerators of the transfer functions 0/6 and 
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7/5 contain a negative linear term yielding a negative static sensitivity 
and causing the l80° initial phase shift. 

For angular frequencies above 10 radians/sec, with the exception of 
u/& (fig- 4 (c)), changing the trim angle of attack has only a slight 
effect on the frequency response obtained, since for this particular 
missile the static stability derivatives have previously been shown to 
be fairly linear. The effect of changing the trim angle of attack would 
therefore be considered negligible in control- system applications where 
the controlling elements are designed to operate at a frequency in excess 
of the high-frequency airframe mode. The reasoning in this case is that 
the overall transient behavior of a control system in combination with 
the airframe would be essentially the same regardless of the trim a, 
since the contribution of the airframe would not change appreciably with 
angle of attack in the region of the control- system operating frequency. 

The curves presented in figure 4 also indicate that the phugoid mode 
is more heavily damped than the high-frequency mode. The values of 
phugoid critical damping ratio given in the first part of table II are 
0.84 at a = 4° and I.96 at a = 10°, compared with values of £ ~ 
for the high-frequency mode. The large value of £ p at a = 10° indi- 
cates that the phugoid mode in this case is nonoscillatory, actually con- 
sisting of two aperiodic modes. This large value of £ p can be attrib- 
uted to the large drag obtained for supersonic flight at this angle of 
attack. The values of undamped natural frequency given in table II where 
the velocity derivatives are included are approximately 0.15 radian/sec 
for the low-frequency mode and 26 radians/ sec for the high-frequency 
mode. The values of the static sensitivity (actual static ratios) of 
the various transfer functions are also given in table II. With the 
exception of a/5, the static sensitivities are infinite at a = 0° and 
finite when a value is assumed for a. In any case the absolute values 
of the 0/6 and 7/ 5 static sensitivities are fairly large (for 
instance, the static value of 9/5 is 29. 1 at a = 4° and -14.8 at 
a = 10° and the infinite value at a = 0° is, of course, due to the 
previously mentioned integration.) These large or negative values of 
static sensitivity, however, may not be too significant since, with the 
assumption of small pert urbat ions, large steady-state values would only 
be applicable for extremely small values of 5. 


Effect of neglecting velocity derivatives .- The frequency responses 
shown in figure 5 are comparable to those shown in figure 4 except that 
in figure 5 the velocity derivatives are neglected (that is, 

= (^ = 0^. The characteristic -equation and static-sensitivity 

values for this case are also summarized in table II. Except for varia- 
tions in static sensitivity, table II, however, does not reveal anything 
very significant concerning the effect of the velocity derivatives. 
Referring again to figures 4 and 5, the variations between the responses 
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shown in figure 5 and the corresponding responses shown in figure 4 are 
mainly in the low-frequency region. A close comparison of the 
0/6 responses shown in figures 4(a) and 5(a) reveals some differences 
in the shapes of the amplitude and phase responses in the low-frequency 
region. The nature of the responses, however, makes the effect of 
including or neglecting the velocity derivatives more readily apparent 
in a comparison of the a/5 responses (figs. 4(b) and 5(b)). In fig- 
ure 5(b) it is shown that changing the angle of attack has very little 
effect on the response. This result arises from the mathematics of the 
transfer function, since the numerator terms practically cancel the 
phugoid quadratic regardless of the angle of attack. However, when the 
velocity derivatives are included, as in figure 4(b), this is not the 
case. Here the amplitude-ratio curves exhibit noticeable change in 
static sensitivity when a value of a is assumed and the phase curves 
show leading characteristics in the low-frequency region. The leading 
characteristics arise from the fact that the numerator terms break at 
lower frequencies than the phugoid mode of motion. 

Effect of Mach number change .- The data presented in table III are 
comparable in every respect to those presented in table II except that 
the values presented are based on a Mach number of 1.2. Variations 
between the corresponding values of the two tables can be noted. Perhaps 
the most significant of these is the value of the phugoid critical damping 
ratio at a = 4° for the case where the velocity derivatives are 
included. In this case f; p = 0.37b and the low value can be attributed 

to a high value of Qd u * Reference to figure 3(n) shows that for the 
a = 4° curve, Cp u « -0.07 or very nearly its negative maximum for this 
angle of attack. The absolute value of Cp u for this case is actually 
greater than C D . (See table I.) Although the curves are not shown, the 
frequency responses for M = 1.2 and a = 4° would exhibit more peaking 
of the amplitude ratio and more rapid phase changing in the region of the 
phugoid undamped natural frequency (0.115 radian/sec) than any of those 
shown for M = 1.6 because of the lower critical damping ratio. 


Three-Degree -of -Freedom Approach to Future 
Automatic-Control Studies 

Since the frequency response obtained with constant flight condi- 
tions is shown to vary with the trim angle of attack, particularly in the 
low-frequency region, the analysis of automatic -control systems which 
operate at low frequencies - for example, altitude or airspeed controls - 
should define the airframe on the basis of three degrees of freedom. The 
transfer-function approach may be essential in preliminary analysis 
employing linear servomechanism theory. Therefore, in analyzing the 
altitude control, for example, depending on the angle-of-attack range 


10 


NACA TN 4192 


anticipated, three or more airframe transfer functions should be com- 
puted for each case with flight conditions and all other system compo- 
nents held constant. Then, depending on the component gain adjustments, 
it is probable that the stability of the system can be adequately pre- 
dicted if no instabilities arise, regardless of the trim angle of attack 
on which the airframe transfer function happens to be based. It is 
obvious that more than one transient solution will be obtained for an 
equilibrium Mach number under these conditions. Although it can only be 
surmised at this time, it is anticipated that the actual transient 
behavior of the entire control system will be within the range of the 
solutions obtained for a number of different trim values of angle of 
attack. In order to obtain the unique solution for a constant set of 
flight conditions, the airframe can be represented by equations of motion, 
wherein some of the coefficients are represented as functions of the vari- 
able a. The airframe defined in this manner is nonlinear but can be com- 
bined with the rest of the automatic-control-system components through 
the use of an analog computer (REAC, for example) to obtain the unique 
solution. This, however, is a lengthy process; nevertheless, it may 
eventually be necessary to determine the validity of the previous 
surmise . 


CONCLUDING REMARKS 


The additional derivatives necessary for solution of the three- 
degree-of-freedom longitudinal equations of motion, namely, the velocity 
derivatives and the slope of the curve of drag coefficient against angle 
of attack, are shown to vary with angle of attack. This variation affects 
the missile frequency responses mainly in the low-frequency region. 
Therefore, in the analysis of a control system which operates at a high 
frequency, the contribution of an airframe with linear static character- 
istics would be essentially the same for all angle-of -attack values. 

A comparison of the frequency responses obtained when the velocity 
derivatives are included with those obtained when the velocity deriva- 
tives are neglected reveals variations in the low-frequency region. The 
angle -of -attack frequency response, for example, exhibits noticeable 
changes in static sensitivity and more pronounced leading phase charac- 
teristics when the velocity derivatives are included than when the 
velocity derivatives are neglected. 

The critical damping ratio was found to be higher for the phugoid 
mode of motion than for the high-frequency airframe mode. For low super- 
sonic Mach numbers, however, where the derivative of the drag coefficient 
with respect to velocity is near its negative peak, the phugoid critical 
damping ratio may assume a fairly low value. 
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For the supersonic canard missile used in this investigation, the 
velocity derivatives reverse sign in the transonic region and have their 
maximum absolute values for supersonic flight in the range of Mach number 
from 1.2 to 1.55* Their absolute values then decrease with further 
increase in Mach number. 

In applying three -degree-of -freedom transfer functions to automatic- 
control- system analysis, more than one airframe transfer function should 
be investigated for an equilibrium Mach number and constant automatic- 
control component gains. Then it is probable that the stability of the 
overall system can be adequately predicted if no instabilities arise, 
regardless of the trim angle of attack on which the airframe transfer 
function happens to be based. It is further surmised that the unique 
transient solution for an equilibrium Mach number will be within the 
range of the solutions based on airframe transfer functions covering a 
range of assumed trim angle-of-attack values. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., February l6, 1954. 
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TABLE I 

C D AM) C L VALUES 


a,, deg 

C D for - 

C L for - 

M = 1.2 

M = 1.6 

M = 1.2 

M = 1.6 

0 

0.05 

0.04 

0 

0 

4 

.063 

.051 

.211 

.182 

10 

.152 

.128 

• 527 

.456 
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TABLE II 


SUMMARY OF TRANSFER-FUNCTION CHARACTERISTIC 


EQUATIONS AND STATIC SENSITIVITIES 


[m 


1.6; sea level; 



Characteristic Equation 


Including velocity derivatives 

a, deg 

, radians/sec 


radians /sec 


0 




26.25 

0.152 

k 

0.132 

0.84 

26.2 

.152 

10 

.207 

I.96 

26.1 

.145 

Neglecting velocity derivatives 

0 

— 



26.25 

.152 

4 

.105 

.83 

26.2 

.152 

10 

.166 

1.31 

26.2 

.152 


Static Sensitivities 


Including velocity derivatives 

a, deg 

0/6 

a/6 

u/5, per radian 

7/6 

0 

CO 

0.9 

00 

00 

4 

29.1 

.5 

-4.2 

28.6 

10 

-14.8 

.5 

-1.7 

-15.3 

Neglecting velocity derivatives 

0 

00 

.9 

00 

00 

4 

41.8 

.9 

-6.7 

4o.9 

10 

-28.9 

.9 

-2.6 

-29.8 
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TABLE III 

SUMMARY OF TRANSFER FUNCTION CHARACTERISTIC 
EQUATIONS AND STATIC SENSITIVITIES 

Jm = 1.2; sea level; x sm = 0.3395] 


Characteristic Equation 


Including velocity derivatives 

cl, deg 

radians/sec 


radians /sec 

i 

0 



22.8 

0.146 

4 

0.115 

0.374 

22.8 

.146 

10 

.182 

1.23 

22.7 

.142 

Neglecting velocity derivatives 

0 



22.8 

.14 

4 

.114 

• 711 

22.8 

.146 

10 

.179 

1.07 

22.7 

.146 


Static Sensitivities 


Including velocity derivatives 

cl, deg 

e/8 

a/5 

u/5, per radian 

y/b 

0 

00 

.8 

00 

00 

4 

8 

.6 

-5.6 

7.4 

10 

-11 

• 5 

-2.2 

-11.5 

Neglecting velocity derivatives 

0 

00 

.8 

00 

00 

4 

26.6 

.8 

-5-7 

25.8 

10 

-16 

.8 

-2.3 

-16.8 


H 

On 



(a) Photograph of model. 

Figure 1.- Supersonic research model. 
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(b) Plan- view sketch of model. 
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Figure 1.- Concluded. 
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Figure 2.- Variation of Cd^ with Mach number for constant values of 

angle of attack. 
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( a ) C D U - 


Figure 3.- Variation of velocity derivatives with Mach number for constant 

values of angle of attack. 
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Figure 3 .- 


Continued . 
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Figure 3*- Concluded 


22 


NACA TN 4192 



(a) 9/6. 


Figure 4.- Three-degree-of- freedom longitudinal frequency responses 

including effect of velocity derivatives for three values of trim a. 
These responses are based on standard sea-level conditions at M = 1.6 
with x sm = 0 . 29 ^'c* 



NACA TN A192 


23 



(t>) a/6. 


Figure k 




Continued 
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(c) -u/8 . 

Figure 4.- Continued 
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Figure 4.- Concluded. 
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(a) 0/5. 

Figure 5 .- Three-degree-of- freedom longitudinal frequency responses 

neglecting effect of velocity derivatives for three values of trim a. 
These responses are based on standard sea-level conditions at M = 1.6 
with x sm = 0 . 29 ^-c. 
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(b) a/5. 

Figure Continued. 
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Figure 5 


Continued 
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(d) 7 / 8 . 


Figure 5*- Concluded, 
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